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Abstract 

We study thermal and chemical equilibration in 'infinite' hadron matter as well 
as in finite size relativistic nucleus-nucleus collisions using a BUU cascade trans- 
port model that contains resonance and string degrees-of-freedom. The 'infinite' 
hadron matter is simulated within a cubic box with periodic boundary conditions. 
The various equilibration times depend on baryon density and energy density and 
are much shorter for particles consisting of light quarks then for particles including 
strangeness. For kaons and antikaons the chemical equilibration time is found to be 
larger than ~ 40 fm/c for all baryon and energy densities considered. The inclusion 
of continuum excitations, i.e. hadron 'strings', leads to a limiting temperature of 
Ts — 150 MeV. We, furthermore, study the expansion of a hadronic fireball after 
equilibration. The slope parameters of the particles after expansion increase with 
their mass; the pions leave the fireball much faster then nucleons and accelerate 
subsequently heavier hadrons by rescattering ('pion wind'). If the system before ex- 
pansion is close to the limiting temperature Tg, the slope parameters for all particles 
after expansion practically do not depend on (initial) energy and baryon density. 
Finally, the equilibration in relativistic nucleus-nucleus collision is considered. Since 
the reaction time here is much shorter than the equilibration time for strangeness, a 
chemical equilibrium of strange particles in heavy-ion collisions is not supported by 
our transport calculations. However, the various particle spectra can approximately 
be described within the blast model. 
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1 Introduction 



Nucleus-nucleus collisions at relativistic and ultrarelativistic energies are studied experi- 
mentally and theoretically to obtain information about the properties of hadrons at high 
density and/or temperature as well as about the phase transition to a new state of matter, 
the quark-gluon plasma (QGP). In the latter deconfined partons are the essential degrees 
of freedom that resolve the underlying structure of hadrons Whereas the early 'big- 
bang' of the universe most likely evolved through steps of kinetic and chemical equilibrium, 
the laboratory 'tiny bangs' proceed through phase-space configurations that initially are 
far from an equilibrium phase and then evolve by fast expansion. These 'specific initial 
conditions' - on the theoretical side - have lead to a rapid development of nonequilibrium 
quantum field theory and nonequilibribrium kinetic theory 0, ^. Presently, semiclassical 
transport models are widely used as approximate solutions to these theories and practi- 
cally are an essential ingredient in the experimental data analysis. For recent reviews we 
refer the reader to Refs. [|, ^ |]. 

On the other hand, many observables from strongly interacting systems are domi- 
nated by many-body phase space such that spectra and abundances look 'thermal'. It is 
thus tempting to characterize the experimental observables by global thermodynamical 
quantities like 'temperature', chemical potentials or entropy ^ |^, |^, |TD|, |TT|. We note, 
that even the use of macroscopic models like hydrodynamics [1^, |T3| employs as basic 
assumption the concept of local thermal and chemical equilibrium. The crucial ques- 
tion, however, how and on what timescales a global thermodynamic equilibrium can be 
achieved, is presently a matter of debate. Thus nonequilibrium approaches have been used 
in the past to address the problem of timescales associated to global or local equilibration 



[14| , p!5| , [Lq , p77| , |Tg, ^ |21|. In view of the increasing 'popularity' of thermodynamic 



analyses a thorough microscopic reanalysis of this questions appears necessary especially 
for nucleus-nucleus collisions at ultrarelativistic energies that aim at 'detecting' a phase 
transition to the QGP. 

In this paper we study equilibration phenomena in 'infinite' hadronic matter using a 
microscopic transport model that contains both hadron resonance and string degrees-of- 
freedom. With this investigation we want to provide insight into the reaction dynamics 
by the use of cascade-like models and also point out some of their limitations. The 
'infinite' hadronic matter is modelled by initializing the system solely by nucleonic degrees 
of freedom through a fixed baryon density and energy density, while confining it to a 
cubic box and imposing periodic boundary conditions during the propagation in time. 
We, furthermore, then study the expansion of the hadronic fireball after equilibration to 
investigate the changes in hadron spectra during the rapid explosion as well as related 
equilibration phenomena in realistic nucleus-nucleus collisions for light and heavy systems. 

Our paper is organized as follows: In Section 2 we briefiy describe the approach 
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employed in our investigations, specify the initial conditions for a finite box with peri- 
odic boundary conditions, present our numerical results and extract various (hadronic) 
equilibration times as well as thermodynamical properties for different initial conditions. 
Section 3 is devoted to the expansion dynamics of the equlibrated fireball and a discussion 
of the related physical phenomena. In Section 4 we analyse reactions of colliding finite 
light and heavy systems and compare our result to a blast model. Section 5 concludes 
our study with a summary. 



2 Equilibration and limiting temperature 



To investigate the equilibration phenomena addressed above we perform microscopic cal- 
culations using the Boltzmann-Uehling-Uhlenbeck (BUU) model of Refs. p3| , |2^ . This 
model is based on the resonance concept of nucleon-nucleon and meson-nucleon interac- 
tions at low invariant energy ^/s [^, adopting all resonance parameters from the Manley 
analysis The high energy collisions - above y/s = 2.6 GeV for baryon-baryon colli- 
sions and y/s = 2.2 GeV for meson-baryon collisions - are described by the LUND string 
fragmentation model FRITIOF ^6\. This aspect is similar to that used in the HSD ap- 
proach ^ |2^ and the UrQMD code . For a detailed description of the underlying 



model at low energy we refer the reader to Ref. [|4 



For later discussions it is essential to realize that the code respects detailed balance 
only for reactions of the type 1 ^ 2 -|- 3 and approximately for l-|-2<->3-|-4Q where 
the numbers 1, ... ,4 are any reaction partners. This implies that in particular at high 
energies, where the string degrees of freedom with their decay to many (> 2) final particles 
becomes important, detailed balance is violated. We will discuss the consequences of this 
violation, which is inherent in all such transport codes, at the appropriate points in the 
following sections. 



2.1 A box with periodic boundary conditions 

In order to study 'infinite' hadronic matter problems we confine the particles in a cubic box 
with periodic boundary conditions for their propagation similar to a recent box calculation 



within the UrQMD model |]TB[. We specify the initial conditions, i.e. baryon density p, 
strange particle density ps and energy density e as follows: first the initial system is fixed 
to Np = 80 protons and Nn = 80 neutrons, which are randomly distributed in a cubic box 
of volume V. The 3-momenta pi of the nucleons in a first step are randomly distributed 
inside a Fermi-sphere of radius pp = 0.26 GeV/c (at po) and in a second step boosted by 

^In the latter case small violations of detailed balance are due to the treatment of t-channel and 
background contributions. 
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i/^cm by a proper Lorentz transformation. Thus the initial baryon density p is fixed as 
p = A/V, A = Np + Nn. The strange particle density is set to zero as in related heavy- ion 
experiments while the energy density is defined as e = E/V, where E is the total energy 
of all nucleons 

i 

The boost velocity jScm is related to the initial energy density e (excluding Fermi motion) 
as 

f^cm = ^1-^ (2) 

using 

£ = IcmPrriN (3) 

with = l/^/l — /3cm- Recall that porrij^ ~ 0.15 GeV/fm^ so that an energy density 
e ~ 1.5 GeV/fm^ at density po corresponds to 7cm — 10, i.e. the SPS energy Tiab — 
185 A-GeV. We thus start with a 'true' nonequilibrium situation in order to mimique 
the initial stage in a relativistic heavy-ion collision. The initial phase represents two 
interpenetrating, (ideally) infinitely extended fiuids of cold nuclear matter. 

We now propagate all particles in the box in the cascade mode (without mean-field 
potentials) using periodic boundary conditions, i.e. particles moving out of the box are 
reinserted at the opposite side with the same momentum. The phase-space distribution 
of particles then can change due to elastic collisions, resonance and string production and 
their decays to mesons and baryons again. We recall that we include all baryon resonances 
up to an invariant mass of 2 GeV and meson resonances up to the </)-meson. According to 
the initial conditions for e and p the factor 7cm in (H) determines if strings are excited in 
the very first collisions. This is the case for 7cm > 1-4 where the early equlibration stages 
are dominated by string formation and decay. 

2.2 Chemical equilibration 

Figure |I] shows the time evolution of the various particle abundances (nucleons A^, A, 
A, TT, ?7, K'^ and mesons) for density p = po (left panel) at different energy densi- 
ties e = 1.1,0.52 and 0.22 GeV/fm^ and for p = 3po (right panel) at e = 3.4,1.57 and 
0.66 GeV/fm^. These initial conditions correspond to bombarding energies Tiab per nu- 
cleon of roughly 100, 20 and 2 A- GeV, respectively. In Fig. |l| (as well as in Figs. |^,§) we 
count all particles which are 'hadronized', i.e. produced by string decay after a formation 
time of Tp = 0.8 fm/c in their rest frame. 
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After several fm/c the number of nucleons decreases due to inelastic collisions that 
produce either baryon resonances or additional mesons. The number of A-resonances 
grows up to a maximum in a few fm/c, since a lot of A's are produced in the first A^A^ 
collisions; their number subsequently decreases with time due to their decay and excitation 
of further resonances or due to reabsorption. The numbers of vr's and r/'s increase very 
fast and reach the equilibrium value within a few fm/c whereas the strange particles 
(AT"*", K", A) require a much longer time for equilibration. 

In Fig. ^we present the time evolution of the particle ratios vr/A^, A/N, A/N, A'^/vr"'", 
K~/tt^, rj/iT for density p = po ai energy densities e = 1.1,0.52 and 0.2 GeV/fm^, while 
Fig. I shows the same particle ratios for density p = 3po at energy densities e = 3.4, 1.57 
and 0.66 GeV/fm^, respectively. The left panels in both plots correspond to the full time 
scale as in Fig. ^ (up to 1000 fm/c), whereas the right panels present in more detail the 
initial phase (up to 30 fm/c). We use the same scale for the |/-axis on the right and left 
panels, so one can easily see that the 7f/N, A/N ratios reach the equilibrium values very 
fast especially at low energy density since the string degrees of freedom here play a minor 
role and pion production basically emerges through A resonance decay. 

The meson-pion ratios {K~^/7t^, K^/tt^, rj/n) at high energies show a decrease in 
the first few fm/c and then an increase again up to the equilibrium values. This is due 
to the fact that the bulk of the strange mesons is produced very early (at high energy 
density) through string formation and decay whereas most of the pions appear later, with 
a delay of several fm/c, as a result of the decay of heavy vector mesons (e.g. p and uj). 
From the right panels of Figs. ^ and ^ one can see that in the initial stage the particle 
ratios containing strange to nonstrange particles - K~^/7r~^, K~ /ti~ , A/N - are still far 
off chemical equilibrium for all energies and densities and the equilibration takes up to a 
few hundred fm/c depending on the energy and baryon density. 

For the higher energies the initial particle production proceeds via the formation and 
decay of string excitations. This leads in particular to a very early onset of strange 
particles (mainly kaons and hyperons) within the first fm/c either due to the initial strings 
or due to secondary or ternary baryon-baryon, meson-baryon and meson-meson induced 
string- like interactions (see the right panels of Figs. ^ and In Ref. |^ it was shown 
that these early secondary and ternary reactions can contribute up to about 50 % of the 
total strange particles obtained in a Pb + Pb reaction at CERN SPS energies and thus 
explain the factor of 2 in the observed relative strangeness enhancement compared to 
p-|-p reactions. This, however, does not imply that chemical equilibrium for the dominant 
strange particles has been achieved in this reaction, as our analysis clearly shows. In the 
later stages, when the system has become, more or less, isotropic in momentum space, 
strange particles can only be further produced by low energy hadronic reactions, which, 
however, have a considerable threshold and are thus strongly suppressed. This explains 
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the long chemical equilibration times for the strange particles first demonstrated by Koch, 
Miiller and Rafelski 0. 

In order to define an overall chemical equilibration time we perform a fit to the particle 
abundances N{t) for pions and kaons as 

iV(t) = iV,Jl-exp(-t/re,)), (4) 

where N^g is the equilibrium limit. The equilibration time Tgg thus corresponds to the 
time t when ~ 63% of Neq is achieved. 

Figure § shows the equilibration time Tgg versus energy density for vr and mesons 
at different baryon densities of l/3po;Po;3po and 6po- We find that the equilibration 
time for pions scales as rj^ ~ 1/p or ~ p, thus we present the curve only for baryon 
density po- Whereas rj^ slowly grows with energy- density, falls steeply with e. This 
marked difference is due to the fact that, on one hand, the kaon production rate increases 
dramatically with ^/s whereas that of the pions, on the other hand, is more fiat. With 
increasing energy thus more strange particles are produced through strings especially from 
the primary collisions with high ^/s and the chemical equilibration is achieved faster. 

In Fig. ^ we have considered an 'ideal' situation, i.e. hadron matter at fixed energy 
and baryon density. In realistic heavy-ion collisions the system goes through the different 
stages due to interactions and expansion. However, as follows from Fig. |^, the equili- 
bration time for strangeness is larger than 40 fm/c for all energy and baryon densities. 
Thus in realistic nucleus-nucleus collisions the chemical equilibration of strange particles 
requires also a time above 40 fm/c which is considerably larger than the actual reaction 
time of a few 10 fm/c or less (cf. Section 4). 

The particle abundances used to extract Teg in Fig. ^ have been calculated without any 
in-medium potentials. In fact, the introduction of attractive potentials (especially for K~) 
will lower the hadronic thresholds and thus increase the scattering rate between strange 
and nonstrange hadrons, whereas the feels some repulsive potential and the trend 
goes in the opposite way. According to our calculations such in-medium modifications (in 
line with Ref. 0) give a correction to the equilibration times by atmost 10 % and 
shortens the K~ equilibration times up to 20 % at density po. 

2.3 Thermal equilibration and limiting temperature 

In this subsection we investigate the approach to thermal equilibration. This is initially 
driven by the very early string phase on the momentum equilibration of the hadronic 
degrees of freedom, when the system is still very far from equilibrium and the energy 
density is sufficiently high. This one can see by looking at the quadrupole moment 
<Q2> =< 2]9^ — pI — Py > of the momentum distribution of all hadrons involved. In 
the left panel of Fig. ^ we present the time evolution of the quadrupole moment < Q2 > 
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for density p = po at energy densities e = 0.22, 0.3, 0.52, 0.8, 1.1 and 1.6 GeV/fni'^. In 
order to take into account the string contributions we have counted here all particles even 
within the formation time. The thin solid lines indicate exponential fits of the form 



<Q2> {t) ^ Aiexp{-t/Tshort) + A2exp{-t/Tiong) 



(5) 



with two equilibration times Tshort and Tiong- 

The right panel of Fig. |^ shows Tshort and Tiong versus energy density e. Whereas 
Tshort — 5 fm/c is roughly independent on e the 'hadronic' equilibration time Tiong in- 
creases with energy density. These results have to be interpreted as follows: in the 
initial nonequilibrium phase the string degrees of freedom are excited and decay accord- 
ing to many-body phase on a short time scale Tshort- The string decays reduce the initial 
quadrupole moment (at high energy density) in time by a significant factor of about 
3 — 4. One can understand the result obtained for Tshort in a rather simple way. Due to 
our prescription of the initialization of the system the first strings on average are formed 
after the time Tcou ~ l/{{p/2)aNN{vNN)) ~ 3 — 4 fm/c for p = po- The strings then 
decay within their formation time Tp ~ 0.8 fm/c giving rise to a significant production of 
transversal momentum. One should point out, that according to these arguments Tshort 
approximately scales like 1/p. Due to Lorentz contraction Tshort is thus considerably 
smaller in a real heavy-ion collision. Hence, string decays provide a very efficient source 
for a strong decrease in longitudinal momentum and production of transverse momentum 
in the very early stage of an ultrarelativistic heavy-ion collision. A decrease (increase) of 
the formation time Tp to 0.5 fm/c (1.5 fm/c) changes Tshort on the scale of 20%, only. 

After string decay, however, the emerging hadronic system still has significantly larger 
longitudinal than transverse momenta - the ratio increases with energy density e - and 
low energy hadronic reactions are less effective in transfering longitudinal to transverse 
momentum or simply in production of mesons. This explains the increase of Tiong with e 
in simple terms. 

From the above analysis it follows that after typical relaxation times of Tshort ~ 5 fm/c 
the momentum unisotropy of hot and dense matter has dropped to e~^ such that one 
might describe the system by simple global thermodynamical variables like temperature 
etc. This thermal equilibrium has to be contrasted with the chemical equilibrium which 
- as we have shown in the preceding subsection - is reached only after much longer times 
(> 40 fm/c for strange particles, for example). 

For the equilibrated system we can extract a temperature T by fitting the particle 
spectra with the Bolzmann distribution 



where Ei = Jpf + mf is the energy of particle i. We note that at the temperatures of 




(6) 
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interest here, the Bose and Fermi distributions are practically identical to a Boltzmann 
distribution. We find that in equilibrium the spectra of all particles can be characterized 
by one single temperature T. This is demonstrated in Fig. ^ where we show the spectra 
of nucleons (A^), pions (vr) and kaons {K'^) as a function of the kinetic energy E — m for 
p = Po at energy densities e = 0.52, 0.8 and 1.6 GeV/fm^ (left panel) and for p = 3po at 
energy densities e = 0.66, 1.57 and 2.85 GeV/fm^ (right panel). Here we have averaged 
the spectra from 950 fm/c to 1000 fm/c in order to decrease the numerical fiuctuations. 
The spectra of N, it, here can be fitted with a single temperature T which increases 
with the energy density e for both baryon densities po and 3po. We note exphcitly that 
the slope of the equilibrium particle spectra does not depend on the formation time Tp. 

In Fig. we display the energy density e versus temperature T for different baryon 
densities p: l/3po (open down triangles), po (full squares), 3po (full dots), 6po (full up tri- 
angles). In order to compare calculations for different baryon densities we have subtracted 
the baryon energy density at rest, i.e. ^ Tn^^p (except for Fermi motion). As seen from 
Fig. the temperature grows with energy density up to a limiting value reminiscent of a 
'Hagedorn' temperature |^0|]. From our detailed investigations we obtain for the limiting 
temperature Tg ~ 150 ± 5 MeV which practically does not depend on baryon density. 
Such a singular behavior of s{T) for T c:^ Tg has also been found in the box calculations 



in Ref. for p = po- Our limiting temperature is slightly higher than that in Ref. flS 



(Tg = 130 ± 10 MeV) due to the different number of degrees of freedom; the model |jT8 



contains more resonances and uses a different threshold for string excitations. Thus, there 
is some phenomenological sensitivity to the hadronic zoo of particles and string thresholds 
employed in the model. 

In Fig. H we show the excitation function for the ratio of string energy density to the 
energy density of the whole system e string I ^ P = Po when the system has equilibrated 
for long times. If the equilibrated system is very dense, lower energy strings are still 
continuously being excited and thus - because of their subsequent decay - the strings 
constitute a stationary portion of the total energy of the system. The relative ratio in the 
energy density increases with e up to a saturation value of ~ 16% and then stays essentially 
constant. This reflects that the system reaches a limiting temperature, since the relative 
amount of string excitations compared to resonance excitations does not change any more, 
whereas the number of strings as well as the number of hadrons produced increases with e. 
This fact one might have guessed since the string production rate in equilibrium depends 
only on the temperature T characterizing the Bose/Fermi distributions in the collision 
terms. In addition, this constant fraction, of course, also intrinsically depends on the 
excitation threshold and on the chosen decay (or formation) time Tp of the strings. 

As pointed out above, the string degrees of freedom play an essential role for particle 
production at high bombarding energies since they describe the continuum excitations 
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of the system. The number of strings created is especially high at the first stages of 
the collision, when the energy of baryon-baryon interactions is close to the initial energy 
\fs. It decreases with time to some constant value which corresponds to the equilibrium 
state. Because of this string-dominance one now has to worry about possible consequences 
of a violation of detailed balance for these degrees of freedom. As already pointed out 
earlier, all hadronic cascade-type approaches use the phenomenological string picture in 
order to describe quantitatively energetic (soft or semi-hard) inelastic reactions above 
some specified A/s-threshold. In such binary hadronic reactions typically many hadronic 
particles and resonances are produced, the number depending on the incident energy 
s/s. The 'back reaction' of these particles produced from decay of an excited string (or 
two strings in the LUND model) leading to the formation of only two energetic hadrons 
again is not considered as it is statistically suppressed and difficult to describe. On the 
other hand, in an 'infinite' matter calculation these back reactions have to be taken into 
account in order to allow for the principle of detailed balance. This is not done here 
as it is technically difficult to handle; it thus represents a potential 'Achilles heal' in a 
thermodynamic analysis. However, for simulating a heavy-ion collision this deficiency is 
not of any major importance since the excitation of strings happens in the first moment 
of the reaction when the phase space is still widely open and no back reaction can occur. 



2.4 Comparison to the statistical model 

In order to investigate the equilibrium behavior of hadron matter we also compare our 
transport (box) calculations with a simple Statistical Model (SM) for an Ideal Hadron Gas 
(IHG) where the system is described by a grand canonical ensemble of non-interacting 
fermions and bosons in equilibrium at temperature T. All baryon and meson species 
considered in the transport model also have been included in the statistical model. 
Our main objective here is to compare our results with the Hagedorn bootstrap picture 
of hadronic matter . 

We recall that in the SM particle multiplicities and energy densities Si are given by 



oo 



{2TxhY I exp [{E, - B.^iB - S,^is)lT] ± 1 ' ^ ^ 





oo 



^ 9i f AnEip'dp 

{2Tihf I exp [{E, - B.fiB - S,fis)/T] ± 1 ' ^ ' 







where Ei = ^p'^ + mf is the energy of particle i, Bi is the baryon charge. Si is the 
strangeness, and is the spin-isospin degeneracy factor. In Eqs. (0),(H) fiB and fis are 
the baryon and strangeness chemical potentials. Here we neglect the electric chemical 
potential (//„ = /Up = /i^) since we consider an isospin symmetric system. Note, however. 
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that in realistic collisions of heavy-ions (like Au + Au) this reduction is no longer fully 
appropriate. For particles with finite spectral width we include in Eqs. (0),(||) the spectral 
functions Pi{m) with the same parametrization for the width as in the transport model, 



oo 



= -T7 / Pi[m]dm I — — ; , 9 







7—1 9 7 

9i f f f A-KEip^dp 

Pi[m)dm / -— — — j . (lUj 



{27Thy J ' y exp [{E, - B^pB - S,ps)/T] ± 1 ' 

The energy density £, baryon density p and strange density of the hole system in 
equilibrium then given as 

^ = ^^iiT,pB,Ps) (11) 

i 

p = ^Bi ni{T, pb, Ps) (12) 

i 

Ps = Yl^i ^iiT,PB,Ps) = 0. (13) 

i 

As 'input' for the SM we use the same e,p and ps as in the box calculations and we 
obtain the thermodynamical parameters - T, ps, ps ^ solving the system of nonlinear 
equations (0),(|D and (|l|). 

Within the SM we find that the temperature increases continuously with energy density 
since the continuum excitations, i.e. the string degrees of freedom, are not included (full 
dots in Fig. P), whereas the box calculation with strings gives the limiting temperature 
(full squares in Fig. ^. Both curves in Fig. ^ have been calculated for density po- 

To reproduce qualitatively our box result within the SM we have to include continuum 
excitations in the statistical model, i.e. a Hagedorn mass spectrum for strings as defined 

by ig 

p^*'-(m) = ^exp(m/TH), (14) 

where Th denotes the 'Hagedorn' temperature. For we use the temperature Tg as 
obtained from the box calculations, i.e. Th = ~ 150 MeV. In ( p!4D p^^^ is a fit parameter 
additionally to T,pB and ps to reproduce e{T) from the box calculations. The string 
multiplicities nf^ are given by 

T °° °° /I 2j 

J J exp [{Ei - BipB - Sips)/T\ ± 1 

where the lowest mass in the string excitation {rrimin) is defined by the string threshold in 
the transport model: m^m = 2.6 — rri]^ GeV for baryon strings and m^m = 2.2 GeV for 
meson strings. In our transport model we include the following strings i: baryon strings 
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B = 1, S = 0, — 1, — 2, — 3, anti-baryon strings B = —1,S = 0,1,2,3, meson strings 
B = 0, S = 0,1 and anti-meson strings B = 0, S = —1. 

Before going over to the actual analysis we point out that the limiting temperature Tg 
from our string model involves somewhat different physics assumptions than the Hagedorn 
model at temperature T^. Tg should not really be identified with the 'Hagedorn' temper- 
ature Th, though close similarities exist. In the Hagedorn picture and for temperatures 
close to Th the abundance of 'normal' hadrons or known resonances stays constant with 
increasing energy density whereas the number and energy density of the (hypothetical) 
bootstrap excitations diverges for T Th- The Hagedorn model thus assumes 'particles' 
of mass m ^ oo to be populated for T Th, that dynamically can be formed in col- 
lisions of high mass hadrons for t ^ oo. In contrast, our string model does not include 
energetic string-string interactions that might produce more massive strings. (There ex- 
ist some phenomenological recipes how to incorporate such interactions [0.) The 'high 
mass' strings decay to hadrons and, because of the detailed balance problem discussed 
in the last subsection, are only repopulated by binary hadron-hadron or hadron-string 
interactions, so that their internal energy is limited and the low-energy hadronic degrees 
of freedom are overpopulated. This leads to the saturation of string-energy to total energy 
(observed in Fig. |^) to a value of ^ 0.16 in contrast to the value of 1 in the Hagedorn 
model. 

This, however, does not imply a fundamental inconsistency for the overall properties 
of the system. In perfect chemical equilibrium, like in the Hagedorn model, more strings 
(or hypothetical resonances) would be excited which, for lower temperatures (e.g. in a 
nearly isentropic expansion of the system like in heavy- ion collisions), would immediately 
decay into a large number of hadronic particles. The violation of detailed balance in our 
case thus physically describes an overpopulation of hadronic particles only in stationary 
equilibrium. The important point, however, is the observation that in either description 
the system at equilibrium can not exceed the critical temperature T^. 

As seen in Fig. R| we achieve agreement of the extended SM and our box calculations 



from Fig. ^ by choosing Th ~ Tg in Eq. (|Tj). In addition, from the extended SM we can 
also define thermodynamical parameters such as the baryon chemical potential /i^. In 
Fig. [ly we present the resulting T — fiB correlation, i.e. temperature T versus baryon 
chemical potential ^b, at fixed baryon densities (in the box calculations) of p = l/3po,Po 
and various energy densities. The open triangles and squares (connected by the dashed 
lines) show the result of the SM without strings at densities l/3po and po, respectively, 
whereas the full triangles and squares (connected by the solid lines) correspond to the 
thermodynamical fit of the box calculations (at l/3po and po) including string excitations. 
The errorbars indicate the uncertainty in the extraction of p^ in the SM; they become 
larger when the system is closer to Th due to the divergence in the energy density integral 
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(|TT|). The arrow at /i^ = indicates the temperature Tg = 150 MeV from our box 
calculations. The full dots in Fig. ^ correspond to chemical freeze-out points extracted 
in a thermodynamical model from hadron abundances the open dots are the thermal 
freeze-out points from the momentum spectra of hadrons and two-particle correlations as 
taken from Ref. |^ . 



Our calculations here are for nuclear matter densities l/3po and po whereas the freeze- 
out points have been extracted from heavy-ion data; the comparison thus can be only 
qualitative. However, one can see the general tendency: if the continuum excitations 
(strings) are not included in the thermodynamical analysis, one can 'extract' much larger 
temperatures at high energy density simply due to the limited number of degrees of 
freedom involved in the model analysis. In this respect our box result is more in line with 
the thermal ('kinetic') freeze-out analysis from Ref. [^] than with the thermodynamical 
analysis from Ref. that is based on particle ratios and thus on chemical freeze-out. 
The point to make is that at higher temperatures, like e.g. the ones obtained for a 
'chemical' freeze-out in Ref. [0, the consideration of continuum excitations does make 
a thermodynamical analysis much less certain than at lower temperatures, like e.g. at 
'thermal' (or kinetic) freeze-out as in Ref. [^, where the continuum excitations do not 
play any significant role. 

In this context we have to mention, furthermore, that a combined experimental anal- 
ysis of particle spectra and HBT radii favors even lower freeze-out temperatures (below 
100 MeV |3^). For these freeze-out conditions the pion density (for fixed charge) 
drops below ~ 10~^ fm'^, i.e. the average distance between two pions (of different charge) 
becomes large than ~ 4.6 fm, which in turn is large compared to their classical interaction 



radius rj = ya^^/vr at all relative momenta between the two pions. Since thermal freeze- 
out temperatures of 90-100 MeV at SPS energies can be considered as a lower bound, the 
'experimental' points in Fig. llO] have to be taken with care. 



3 Expanding hadronic fireballs 

In realistic nucleus-nucleus collisions the system rapidly expands after the possible forma- 
tion of a hot hadronic fireball. The final hadronic spectra can be changed substantially 
during this expansion phase, i.e. the temperature extracted from the experimentally 
observed slopes of the spectra also contains information about the nuclear expansion 
dynamics. 

To investigate the expansion of the hadronic fireball we initialize the system in a box 
with periodic boundary conditions - as described above - and propagate the system up 
to 500 fm/c, when equilibrium is reached. Afterwards we let the system expand without 
boundary conditions. Even though this is an idealized description of the expansion phase 
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during a heavy-ion collision we hope to learn from this scenario how the expansion stage 
changes the picture of perfect thermal equilibrium (for an analysis of an actual collision 
see the discussion in the following section). 

In Fig. |ll| we present the time evolution of the various particle abundances (nucleons 
N, A, A, TT, and mesons) during the expansion for density p = po (left panel) 
at different energy densities e = 0.22,0.3 and 1.1 GeV/fm'^ and for density p = l/3po at 
e = 0.84 GeV/fm^ (upper part in the right panel), for p = po at £ = 1.6 GeV/fm^ (middle 
part in the right panel) and for p = 3po at e = 3.4 GeV/fm^ (lower part in the right 
panel). The number of stable particles [N, A,tt, , K~) increases during the expansion 
up to some asymptotic value due to string and heavy resonance decay as well as inelastic 
interactions. One can see that the asymptotic values are reached after a few 10 fm/c from 
the beginning of the expansion (depending on the initial energies and baryon densities) 
which is comparable to the actual reaction time in heavy-ion collisions (cf. Section 4). 

In Fig. |l^ we show the spectra of nucleons (A^), pions (vr) and kaons {K'^) versus the 
kinetic energy E — m for p = po at energy densities e = 0.22, 0.3 and 1.1 GeV/fm^ before 
the expansion - averaged over time from 450 fm/c to 500 fm/c - (left panel) and after 
the expansion - averaged from 580 fm/c to 600 fm/c - (right panel). For completeness 
in Fig. |1^ we present the result for p = l/3po at e = 0.84 GeV/fm^ (upper part), for 
p = Po at £ = 1.6 GeV/fm^ (middle part) and for p = 3po at e = 3.4 GeV/fm^ (lower 
part). In the left panels the systems are in equilibrium; the N,7i,K^ spectra show a 
common temperature T whereas after the expansion the slopes of the particle spectra are 
different; the nucleon spectra are much harder than the pion spectra, i.e. the apparent 
temperature of particles (after the expansion) increases with the mass m. This effect 
is illustrated in Fig. 0, where we show the apparent slope T versus m for n,K~^,N for 
p = Po at different energy densities: e = 0.2 GeV/fm^ (full up triangles), e = 0.3 GeV/fm^ 
(full squares), e = 1.1 GeV/fm^ (full dots), e = 1.6 GeV/fm^ (full diamonds); and for 
p = 3po at £ = 3.4 GeV/fm'^ (open down triangles). The arrow indicates the limiting 
temperature T, = 150 MeV before the expansion. One can see from Figs. p!^-[l^ that the 
'expansion' temperature of particles increases also with the energy density. However, if 
during the equilibration phase the system reaches T,, the 'expansion' temperatures for 
different particles show a universal behaviour, i.e. practically do not depend on the energy 
e as well as on the baryon density p. This phenomenon is due to the fact that close to 
Tg the initial hadron velocity distributions, reflected in the particle momentum profile, 
become similar for all e and p in equilibrium. 

In order to investigate the origin of the enhancement in the particle slope during the 
expansion we have performed several illustrative calculations: at 500 fm/c - after the 
system has achieved equilibrium - we i) let all resonances and strings decay; in this case 
we find that the slopes do not change as compared to the equilibrium phase, ii) we let 
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the system expand without interactions (allowing only decays) and find that the slopes 
slightly decrease in comparison to the equilibrium phase. Both examples indicate that 
the slope enhancement stems basically from multiple interactions of the particles in the 
initial stages of the expansion phase. 

For analyzing the expansion fiow phenomena we have performed a fit of the particle 
spectra (after expansion) using the blast model of Siemens and Rasmussen [0. In this 
model all particle spectra are described by a universal formula with common thermal 
freeze-out parameters, i.e. a temperature T of the fireball and a radial- flow velocity (3: 



Ai exp ( — ^ ) _^ _ Z. cosh ai 



(16) 



where 7 = (1 - /5^)"^/^ a = -fPpi/T. Here Ei ,Pi are the total energy and momentum of 
the considered particle i while are normalization factors. 

We now try to describe the final particle spectra after the expansion by Eq. (0) with 
common freeze-out parameters T and f3. In Fig. ^ we show the result of our least-squares 
fit, using the MINUIT method |3^, for the energy densities e = 0.3 GeV/fm^ (upper part) 
and 1.1 GeV/fm'^ (lower part) at p = po- The left panel shows the contour plots for the 
parameter errors in the T — [3 plane (for the X^umai + 1 level); the dot-dashed lines stand 
for nucleons (A^), the solid lines for pions (vr) and the dashed lines for kaons {K'^). The 
full symbols indicate the 'best' values for T and (3 according to the criteria (squares for 
A^, dots for vr and triangles for K^). The thin solid lines in the right panel demonstrate 
the fit of the particle spectra within the optimal parameters from MINUIT. 

Since the particle spectra cover several orders of magnitude and the low energy points 
contribute to with a larger weight than those at high energy, we use the logarithmic 
method to give a higher weight to the tail of the spectra in the fitting procedure, i.e. we 
minimize Xki = Z](ln/(xj) — ln/o(xi))^, where /o represent the 'experimental' data (i.e. 

i 

the results of our box calculations), / is the value of the fit ([T6|) at point Xj. 

One can see from Fig. that the 'best' parameters T and (3 (as well as the contours for 
the parameter errors) are quite different for N , ir and especially for e = 0.3 GeV/fm^. 
So we do not find (within the 'optimal' x^) common freeze-out parameters for all spectra 
simultaneously. This is similar to an analysis of experimental spectra by Peitzmann et 
al. |3^. For all particles we obtain different values for (3 and much lower temperatures T 
than that of the initial fireball: Tin = 107 MeV for e = 0.3 GeV/fm^ and Ti„ = 145 MeV 
for e = 1.1 GeV/fm^. 

Since especially the pion spectra contain large contributions from resonance decays 
at low E ~ m, we have also performed fits excluding the pion spectra for E — m < 0.4 
GeV. This procedure essentially gives lower f3 parameters and higher values for T (open 
circle in the upper panel). However, the low energy cut-off is an additional parameter 
that allows to 'extend' the (/5, T) values to a wider range. 
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Thus our analysis indicates that the final particle spectra do not allow a reliable recon- 
struction of freeze-out parameters within the collective flow model (16). The parameters 
T and (3 obtained from the fit are very sensitive to the low energy shape of the hadron 
spectra or low energy cut-off applied since this region contributes with the largest weight 
to the minimization. On the other hand, a global 'eye' fit with the parameters given 
by the 'star' for all hadrons considered gives a quite reasonable overall description of the 



spectra (dashed lines, r.h.s of Fig. [TSD, too. A very accurate deduction of one single overall 
fit for all hadrons by a common temperature (T) and flow velocity {(3) parameter (see, 
e.g., 1^ and references therein) seems to us thus rather ambiguous. In particular, such an 
analysis may indicate that thermal equilibrium has been reached to a much larger extent 
than is actually true. 

The particle flow effect due to the expansion is demonstrated in Fig. ^ where we show 
the velocity distributions dN/djS for nucleons (A^), pions (vr) and kaons {K~^) for p = l/3po 
at £ = 0.84 GeV/fm^ (upper part), for p = po at e = 1.1 GeV/fm^ (middle part) and for 
p = 3po at £ = 3.4 GeV/fm'^ (lower part). The left panel shows the dN/d(3 distribution at 
equilibrium whereas the right panel corresponds to dN/d(3 after the expansion phase. On 
can see that the average velocity of the particles decreases with the mass; the pions are 
much faster than the nucleons. They thus leave the reaction zone at the initial stage of the 
ongoing and rapidly evolving expansion with a higher velocity and accelerate the slower 
hadrons that 'feel' the 'pion wind' by the multiple interactions We recall that the 
pion density is very high especially at high e such that practically all other hadrons are 
shifted in direction of large (3. The same effect is shown in Figs. pr2| , p!3| by the enhancement 
of the slope of the nucleon spectra due to the expansion. 



4 Reactions of finite systems 

In this Section we turn to realistic nucleus-nucleus collisions with the BUU transport 
model. We have learned from our analysis in the previous Section that even by starting 
from an idealized scenario of perfect thermal equilibrium, a rapid expansion stage makes 
the extraction of one global temperature T and one global expansion parameter (5 quite 
ambiguous. We thus expect this to become even worse for the true situation of a rela- 
tivistic heavy-ion collision, where a perfect equilibrium state at some intermediate stage 
cannot really be assumed. Also, we note that even for a very heavy system like Pb+Pb 
the fraction of effective surface layer (~ AttRIj^X) to total volume is still quite sizeable, 
resulting in a continuous emission or evaporation of particles from the outer layers be- 
fore a global freeze-out of bulk particles occurs. (For Pb-|-Pb reactions at SPS energies - 
combining a hydrodynamical evolution with a nonequilibrium picture of surface emission 
- it has indeed been shown that at least 25% of all particles are continuously evaporated 
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before a global freeze-out has occurred 



In Fig. we show the time evolution of the particle abundances (nucleons A^, A, 
A, TT, rj, K^) for central collision of the light system ^^C + ^^C (upper part) and heavy 
system ^^''Au + ^^"^Au and ^osp]^ _|_ 208p]^ (lower part) at the low energy of 1 A-GeV 
(left panel) and the high energy of 100 and 160 A-GeV (right panel). The number of 
nucleons decreases in a few fm/c due to the inelastic collisions, whereas the number of 
A-resonances increases accordingly. At low energy the 77-mesons and strange particles (we 
disregard strange particles for C + C at 1 A-GeV due to the low statistics) appear with 
a delay of a few fm/c due to the fact that they are basically produced from resonance 
decays (the same as pions) or from secondary pion-baryon collisions, whereas at high 
energy they appear earlier due to the primary production mechanism through the string 
formation and decay. 

As seen from Fig |1^ the reaction time Treac for 1 A-GeV is ~ 20 — 30 fm/c, whereas 
for high energies Treac is shorter due to a faster expansion - r^eac — 10 — 20 fm/c. It has 
been shown in Section 2.2 that the chemical equilibration of hadronic matter under 'ideal' 
conditions (box without expansion) requires a quite long time, e.g. the equilibration time 
Teq for strange particles has been found to be larger than 40 fm/c for all energies and 
densities (cf. Fig. ^). In realistic central nucleus- nucleus collisions, such as Au + Au, the 
system expands rapidly (depending on the energy) after the compression and formation 
of the hadronic fireball. The number of interactions, which is the dynamical origin for 
equilibration, decreases correspondingly very fast with time; after a few 10 fm/c the 
particles are moving practically freely. Thus, the reaction time even for central Au + Au 
collisions is much shorter than the time required for strangeness equilibration: Treac Teq. 

The thermal equilibration time in this energy range around 0.25 GeV/fm^, as obtained 
from the box calculation, is about 5-7 fm/c (see Fig. Notice, however, that this 
calculation - because of its periodic boundary conditions - probably underestimates the 
equilibration time. Indeed, studies of the longitudinal and transverse temperatures {Tl 
and Tt, resp.) |jl6|, ^ have shown that full thermal equilibrium is reached only in the very 
late expansion phase, when the density p has dropped already below its saturation value. 
After a period of 10 fm/c (after first contact) one still finds ~ 1.5 Tt, i.e. an anisotropy 
of about 40 %, considerably more than indicated in Fig. |^. At the bombarding energy of 
160 A-GeV we find a rapid decrease of the quadrupole moment in the momentum space of 
all hadrons by about a factor of 3 at the scale of 5 fm/c leading to longitudinally expanding 
matter. In view of Fig. |^ this stretched ellipsoid in momentum space becomes isotropic 
only on the scale of 10-20 fm/c since low energy hadronic reactions are less effective for 
equilibration. Since this 'hadronic' equilibration time is larger than the reaction time 
for Pb + Pb at 160 A-GeV a substantial anisotropy remains in the hadron momentum 
distributions after the collision. 
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Contrary to the box case the d?N/ dp^ spectra for reahstic nucleus-nucleus collisions do 
not follow the simple exponential behaviour (P) due to the strong longitudinal expansion; 
especially at high bombarding energy the particle spectra show the specific 'banana' shape 
(reflecting the pp spectra at high energies). In order to exclude this simple dynamical effect 
related to the longitudinal expansion, we present in Fig. ^ (right panel) the transverse 
mass spectra l/m^dN/drriT at mid-rapidity (—0.5 < ycm < 0.5) versus mT — m for central 
Au -|- Au collisions at 1 A-GeV (upper part) and for central Pb + Pb at 160 A-GeV 
(lower part) calculated at the end of the reaction. The my-spectra show an exponential 
behaviour (excluding small itlt), however, with different slopes which can not be 
associated directly with a temperature of a hot fireball formed at the intermediate stages 
of the reaction. 

In line with Section 3 we have performed a fit within the blast model (|T^ within 
an interval of unit rapidity around midrapidity using MINUIT. The results of the fit are 
displayed in Fig. |18[ The full symbols (squares for N, dots for vr and triangles for K^) 
correspond to the 'best' values for T and (3 according to the criteria. The thin solid lines 
in the right panel demonstrate the fit of the rriT spectra within the optimal fit parameters 
(we obtain a smaller within the linear method (x^ = J2{f{xi) — foixi))"^), which 

i 

provides a better description of the low niT spectra). 

Similar to Section 3 (cf. Fig. ITH) we obtain (within the 'optimal' criterium) quite 
different freeze-out parameters for A^, vr and spectra. In order to exclude the influence 
of A- (and other resonance) decays on the pion spectra and to investigate the sensitivity of 
the freeze-out parameters to the low energy cuts applied, we performed a fit of the particle 
spectra using the following cut-offs: rriT — m > 0.2 GeV (open symbols) and itlt — m > 
0.4 GeV (open symbols with crosses inside) for Au -|- Au at 1 A-GeV; rriT — m> 0.4 GeV 
(open symbols) and mT — m > 0.5 GeV (open symbols with crosses inside) for Pb + Pb at 
160 A-GeV. As seen from the left panel of Fig. [18| the implementation of the low my cut- 
off leads to a substantial shift of the 'optimal' MINUIT parameters (3 and T especially 
for pions. The (3, T values for the different spectra move towards to each other when 
discarding the low points. For Au -|- Au at 1 A-GeV our f3 and T parameters agree 



with those extracted by the TAPS collaboration using the blast model (star in the 
upper left plot). Here we have to mention that the cut-off — m > 0.4 GeV has been 
applied in the experimental analysis, too |^T|. For the Pb -|- Pb spectra at 160 A-GeV our 
freeze-out parameters are similar to those from Kampfer |42[ (T = 120 MeV, (3 = 0.43; 
star in the lower left plot). The dashed lines in the right panel of Fig. show the fit 
to the particle mx spectra for the (3, T values corresponding to the 'stars' from the left 
panel. Again this 'eye' fit gives a reasonable description of the spectra (except of the very 
low mx part). 

Here we have to mention again that the extraction of freeze-out parameters from the 
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experimental data is very sensitive to the details of the thermodynamical model applied 
as well as to the observables considered. For example, the analysis of SIS data at 1.0 
GeV from Ref. |TT] gives thermal freeze-out parameters - T ~ 52 MeV and /3 ~ 0.4. At 
SPS energies the chemical freeze-out temperature extracted in Ref. |^ from the thermal- 
analysis of particle ratios is T ~ 168 MeV, whereas the analysis of particle spectra and 
two-particle correlations (HBT data) ^3], Q provides a much lower thermal freeze-out 
temperature T ~ 90 — 95 MeV. For a survey different freeze-out parameters the reader is 
referred to Fig. 4 of Ref. fjl ). 

In view of the various uncertainties inherent in the extraction of the thermal freeze-out 
parameters we conclude that a full, i.e. thermal and chemical, thermodynamical equilib- 
rium at freeze-out cannot be deduced from such an analysis. 



5 Summary 

In this paper we have performed a systematic study of equilibration phenomena and 
equilibrium properties of 'infinite' hadronic matter as well as of relativistic nucleus-nucleus 
collisions using a BUU transport model that contains resonance and string degrees-of- 
freedom. The 'infinite' hadron matter is modelled by initializing the system at fixed 
baryon density, strange density and energy density by confining it in a cubic box with 
periodic boundary conditions. 

We have shown that the equilibration times Tgg for different particles depend on baryon 
density and energy density. The time Teq for non-strange particles is much shorter than 
for particles including strangeness; for kaons and antikaons the equilibration time is found 
to be larger than ~ 40 fm/c for all baryon and energy densities considered. The overall 
abundance of the dominant strange particles (kaons and A's) being produced and obtained 
within the BUU cascade model for heavy-ion collisions can therefore not be described by 
assuming a perfect chemical equilibrium as strangeness is typically still undersaturated to 
a quite large extent. We mention taht transport model calculations like ours can describe 
the yield and spectra of the produced nonstrange hadrons as well as , K~, A yields quite 
well at SPS energies On the other hand, at AGS energies the measured K'^/n^ 



ratio in central Au + Au collisions is underestimated by about 30% However, we 

have to point out that the more exotic strange particles (like the measured antihyperon 
yields of Ref. [Q) can by far not be explained within such standard hadronic multiple 
channel reactions. These hadronic data possibly point towards new physics. 

We have, furthermore, shown that thermal equilibrium is established quickly, within 
about 5 fm/c at SIS energies and samewhat larger times at high energies. The inclusion 
of continuum excitations, i.e. hadron 'strings', leads to a limiting temperature of ~ 
150 MeV in our transport approach which practically does not depend on the baryon 
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density and energy. We have compared our results with the statistical model (SM), which 
contains the same degrees of freedom and the same spectral functions of particles as our 
transport model. We found that the limiting temperature behaviour can be reproduced 
in the statistical model only after including continuum excitations of the Hagedorn type, 
otherwise the fireball temperature extracted from the particle abundances and spectra is 
overestimated substantially. 

Close to the critical temperature T^, the hadronic energy densities can increase to a 
couple of GeV/fm^. From lattice QCD calculations one expects that a phase transition to 
a potentially deconfined QGP state should occur. Referring to the limiting temperature 
Ts ~ 150 MeV obtained, a QGP should be revealed and clearly distinguished from a 
hadronic state of matter if one can unambiguously prove the existence of an equilibrated 
and thermal phase of strongly interacting matter with temperatures exceeding, e.g., 200 
MeV. The best candidates are electromagnetic probes, either direct photons or dileptons. 
On the other hand these are also 'contaminated' by hadronic background and/or pree- 
quilibrium physics. So far no thermal electromagnetic source with temperatures larger or 
equal than 200 MeV has been clearly identified. 

We have also studied the expansion of the equilibrated hadronic fireball and found that 
the slope parameters of the particles after expansion increase with their mass; the pions 
leave the fireball much faster than nucleons and accelerate heavier hadrons by rescattering 
('pion wind'). If the system before expansion is close to the limiting temperature T^, the 
slope parameters for all particles after expansion practically do not depend on energy 
and baryon density. This is due to the fact that the particle velocity distributions in 
equilibrium do not change any more for T ^ Tg. We have fitted the resulting spectra 
within the blast model of Siemens and Rasmussen. Our analysis shows a strong sensitivity 
of the (/?, T) parameters on the spectral shape at low energy (or a low energy cut-off) 
so that no reliable parameter determination can be reported. However, a global 'eye' fit 
with 'average' T) parameters describes the data reasonably well. 

Additionally, we have considered the equilibration in realistic nucleus-nucleus collisions 
of light (C + C) and heavy (Au + Au and Pb + Pb) systems. The T) parameters 
extracted from our calculations for Au + Au at 1 A-GeV agree with those extracted from 



the TAPS collaboration 41 1 and for Pb + Pb at 160 A-GeV with the parameters from 



Ref. Here the reaction time is a few 10 fm/c and decreases with the initial energy 
due to the fast expansion. Since the reaction time is much shorter than the equilibration 
time for strangeness, a chemical equilibrium of strange particles in heavy-ion collisions is 
not supported by our transport calculations. Although again simple fits within the blast 
model provide a decent parametrization of our transport results for the differential particle 
spectra (Fig. 0) a deduction of global parameters for thermal freeze-out is again found to 
be rather ambiguous, especially when considering also the lower momentum contributions 
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of the various particle spectra. 
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Figure 1: Time evolution of the various particle abundances (nucleons A^, A, A, vr, rj, 
and K~ mesons) for density p = Po (left panel) at different energy densities e = 1.1,0.52 
and 0.22 GeV/fm^ and for p = 3po (right panel) at e = 3.4, 1.57 and 0.66 GeV/fm^. 
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Figure 2: Time evolution of particle ratios vr/A^, A/N, A/N, K^/n^, K^/n^, r^/vr for 
density p = po at energy densities e = 1.1,0.52 and 0.2 GeV/fm'^. The left panel shows 
the time scale up to 1000 fm/c, whereas the right panel demonstrates the initial stage up 
to 30 fm/c. 
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Figure 3: Time evolution of particle ratios ir/N, A/N, A/N, K^/n^, K^/n^, r^/vr for 
density p = 3po at energy densities e = 3.4, 1.57 and 0.66 GeV/fm^. The left panel shows 
the time scale up to 1000 fm/c, whereas the right panel demonstrates the initial stage up 
to 30 fm/c. 
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Figure 4: Equilibration time Teq versus energy density e for tt and mesons at different 
baryon densities l/Spo^Po^Spo and 6po- 
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Figure 5: Left panel: time evolution of the quadrupole moment < Q2 > for density 
p = Po at energy densities e = 0.22,0.3,0.52,0.8,1.1 and 1.6 GeV/fm^. The thin sohd 
lines indicate the fit of < Q2 > by two exponentials (|^). The parameters Tghort and Tiong 
are shown in the right panel as a function of the energy density e. 
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Figure 6: The spectra of nucleons (A^), pions (tt) and kaons {K~^) as a function of the 
kinetic energy E — m for p = po at energy densities e = 0.52,0.8 and 1.6 GeV/fm^ (left 
panel) and for p = 3po at energy densities e = 0.66, 1.57 and 2.85 GeV/fm^ (right panel). 
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Figure 7: The energy density e — rriNp versus equilibrium temperature T for different 
baryon densities p: l/3po (open down triangles), po (full squares), 3po (full dots), 6po (full 
up triangles). 
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Figure 8: The excitation function for the ratio of string energy density to the energy 
density of the hole system e string/^ at p = po- 
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Figure 9: The energy density versus equilibrium temperature T for baryon density p = Pq. 
The full dots correspond to the statistical model (SM) without strings, the full squares 
show our box calculations including string degrees of freedom, while the solid line shows 
the result from the extended SM including a Hagedorn mass spectrum for strings. 
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Figure 10: The T — fiB phase correlation, i.e. temperature T versus baryon chemical 
potential fiB- The open triangles and squares (connected by the dashed lines) show the 
result of the statistical model without strings (standard SM) fitted to out box calculations 
at densities l/3po and po, respectively, whereas the full triangles and squares (connected 
by the solid lines) correspond to the thermodynamical fit of the box calculations (at l/3po 
and Po) including string excitations (extended SM). The arrow at p^ = indicates the 
limiting temperature Tg = 150 MeV from our box calculations. The full dots correspond to 
the chemical freeze-out points from Ref. while the open dots are the thermal freeze-out 



points from Ref. ||32|| . 
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Figure 11: Time evolution of the various particle abundances (nucleons A^, A, A, vr, 
and K~ mesons) during the expansion (starting at t = 500 fm/c) for density p = po 
(left panel) at different energy densities e = 0.22,0.3 and 1.1 GeV/fm^ and for density 
p = l/3po at e = 0.84 GeV/fm^ (upper part in the right panel), for p = po at e = 1.6 
GeV/fm^ (middle part in the right panel) and for p = 3po at e = 3.4 GeV/fm^ (lower 
part in the right panel). 
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Figure 12: The spectra of nucleons (A^), pions (vr) and kaons (K^) versus the kinetic 
energy E — m for p = po at e = 0.22,0.3 and 1.1 GeV/fni'^ before the expansion (left 
panel) and after the expansion (right panel). 
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Figure 13: The spectra of nucleons (A^), pions (tt) and kaons {K~^) versus the kinetic 
energy E — m before (left panel) and after (right panel) expansion for p = l/3po at 
e = 0.84 GeV/fm^ (upper part), for p = po a.t e = 1.6 GeV/fm^ (middel part) and for 
p = 3po at e = 3.4 GeV/fm^ (lower part). 
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Figure 14: The spectral slope T after expansion versus the hadron mass m for vr, , N 
at p = po and different energy densities: e = 0.2 GeV/fm^ (up full triangles), e = 
0.3 GeV/fm^ (full squares), e = 1.1 GeV/fm^ (full dots), e = 1.6 GeV/fm^ (full diamonds); 
for p = 3po at £ = 3.4 GeV/fm'^ (down open triangles). The arrow indicates the limiting 
temperature Tg ~ 150 MeV before the expansion. 
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Figure 15: Left panel: the contour plots for the parameter errors in the T — /3 plane; dot- 
dashed lines: nucleons (A^), solid lines: pions (vr) and dashed lines: kaons (K^) for the 
energy densities e = 0.3 (upper part) and 1.1 GeV/fm^ (lower part) at p = po- The full 
symbols indicate the 'optimal' parameters T and (3 (squares for A^, dots for vr and triangles 
for K~^). The open dot (upper left plot) reflects /?, T for the pion spectra including the 
cut-off E — m > 0.4 GeV. Right panel: the full symbols (squares for N, dots for vr and 
triangles for K~^) are the box calculations (for the same p and e as in the left panel). The 
thin solid lines show the fit of the particle spectra with the 'optimal' T and f3 parameters 
from the left panel; the dashed lines correspond to the 'eye' fit with the average f3 and T 
parameters given by the stars from the left panel. 
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Figure 16: The velocity distributions dN/d[3 for nucleons (A^), pions (tt) and kaons {K^) 
for p = l/3po at e = 0.84 GeV/fm^ (upper part), for p = pq a.t e = 1.1 GeV/fm^ (middle 
part) and for p = 3po at e = 3.4 GeV/fni'^ (lower part). The left panel shows dN/d(3 at 
equilibrium whereas the right panel corresponds to dN / df3 after the expansion phase. 
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Figure 17: Time evolution of the particle abundances (nucleons A^, A, A, tt, t], K 
central C + C collisions (upper part) and Au + Au and Pb + Pb collisions (lower 
at 1 A-GeV (left panel) and 160 A-GeV (right panel). 
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Figure 18: Left panel: tlie full symbols indicate the 'optimal' parameters T and P (squares 
for A^, dots for vr and triangles for K~^) obtained by exploring the full rriT spectra in the 
fitting procedure. The open symbols correspond to (3, T for the cut-off rriT — m > 0.2 GeV 
(upper part) and 0.4 GeV (lower part); the open symbols with crosses inside indicate the 
/?, T parameters for the cut-off niT — m > 0.4 GeV (upper part) and 0.5 GeV (lower part). 
The stars corresponds to the (3, T parameters from Ref. for Au + Au at 1 A- GeV and 
from Ref. |^ for Pb + Pb at 160 A- GeV. Right panel: Transverse mass (m-r) spectra of 
nucleons (A^), pions (vr) and kaons {K~^) for central Au + Au collisions at 1 A-GeV (upper 
part) and central Pb + Pb collisions at 160 A-GeV (lower part) for —0.5 < i/cm ^ 0.5: 
the full symbols are the transport calculation; the thin solid lines show the spectra 
for the 'optimal' T and (3 parameters from the left panel; the dashed lines show the fit 
with the (3, T values corresponding to the 'stars' from the left panel. 



40 



